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We study spin j hadronic states in quenched lattice QCD to search for a possible S = +1 pen- 
taquark resonance. Our work is the first systematic lattice QCD study which properly carries out 
the following analyses: (1) the careful extraction of the first two low energy states with very high 
statistics and the variational method and (2) the study of volume dependences of eigenenergies 
and spectral weights to distinguish resonance states from scattering states. 
Simulations are carried out on 8 3 x 24, 10 3 x 24, 12 3 x 24 and 16 3 x 24 lattices at j3=5.7 with the 
standard plaquette gauge action and the Wilson quark action. Our result indicates the existence of 
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1. Introduction 

After the first discovery ^ of © + (1540), identifying the properties of the particle is one of 
the central problems in hadron physics. While the isospin of + is likely to be zero, the spin and 
the parity and the origin of its tiny width still remain open questions. In spite of many theoretical 
and experimental studies on & + [Q], the nature of this exotic particle, including the very existence 
of the particle, is still controversial. In such a situation, the lattice QCD calculation is considered 
as one of the most reliable methods, which is an ab initio method directly based on QCD. 

Up to now, several lattice QCD studies have been reported, which mainly look for pentaquarks 
in various different ways [§, §]. However, the conclusion has not been unfortunately settled yet. 
The difficulties in the spectroscopy calculation with lattice QCD generally arise from systematic 
errors due to the discretization, the chiral extrapolation, the quenching effect, the finite volume 
effect and the contaminations from higher excited-states. The difficulty specific to the current 
problem is that the signal of @ + is embedded in the discrete spectrum of NK scattering states 
in finite volume. In order to verify the existence of a resonance state, one needs to isolate the 
first few low energy states including the lowest NK scattering state, identify a resonance state and 
study its volume dependence which can distinguish itself from other scattering states. Ideally, 
one should extract multistates from a high statistics unquenched calculation for several different 
physical volumes, where both the continuum and the chiral limits are taken. However, there are no 
lattice QCD study which performs all these steps so far, due to the enormous computational costs 
needed. We therefore concentrate ourselves on the analyses using rather heavy quarks on coarse 
quenched lattices but with a good statistics and a proper separation of states. 

In this report, we study (/,/) = (0, j) channel in quenched lattice QCD to search for possible 
resonance states. We prepare thousands of gauge configurations in order to single out the signals 
of a resonance state with good statistics. We adopt two independent operators with 1 = and 
J = I and diagonalize the 2x2 correlation matrices for all the combinations of lattice sizes and 
quark masses to extract the 2nd-lowest state slightly above the NK threshold in this channel. After 
the careful separation of the states, we investigate the volume dependence of the energy as well 
as the spectral weight ^ of each state so that we can distinguish the resonance signal from the 
background signals of NK scattering states. Our calculation (For the details, please see our recent 
paper, Ref. [§].) gives the first systematic study which performs 

• the analysis using correlation matrices with very high statistics, which enables us to isolate 
the first few low energy states including the lowest NK scattering state 

• the study of the volume dependences of eigenenergies and spectral weight factors with sev- 
eral different physical volumes, which enables us to distinguish a resonance state from scat- 
tering states 

in a proper way. 

2. Setup 

We adopt the variational method using correlation matrices constructed from independent op- 
erators [Q] in order to isolate the signal of a resonance state from the NK scattering states. The sim- 
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ulations are carried out on four different sizes of lattices, 8 3 x 24, 10 3 x 24, 12 3 x 24 and 16 3 x 24 
with 2900, 2900, 1950 and 950 gauge configurations, respectively, using the standard plaquette 
(Wilson) gauge action at /3 = 5.7 and the Wilson quark action. The hopping parameters for the 
quarks are (k^, JC,)=(0.1600,0.1650), (0.1625,0.1650), (0.1650,0.1650), (0.1600,0.1600) and 
(0.1650,0.1600), which correspond to the current quark masses (m U}C i,m s ) ~ (240, 100), (170, 100), 
(100, 100), (240,240) and (100,240), respectively in the unit of MeV. The lattice spacing a de- 
termined from the Sommer scale is about 0.17 fm, which implies the physical lattice sizes are 
1.4 3 x 4.0 fm 4 , 1.7 3 x 4.0 fm 4 , 2.0 3 x 4.0 fm 4 and 2.7 3 x 4.0 fm 4 . For the con-elation matrices, 
we adopt two wall-type operators © wa u(?) and © wa u(0 1@L whose spinor and color structures are 
the same as those adopted in the first paper by Csikor et al. Periodic boundary conditions are 
taken in all directions for the gauge field, whereas we impose periodic boundary conditions on the 
spatial directions and the Dirichlet boundary condition on the temporal direction for the quark field 
in order to avoid possible contaminations from those propagating beyond the boundary at t = in 
(anti)periodic boundary conditions [Q], which is peculiar to pentaquarks and has not been dwelled. 
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Figure 1: The black (gray) filled-squares denote the lattice QCD data of the 2nd-lowest state in (I, J ) = 
(0, j ) channel extracted with A^ t — 1 (A^ t ) data plotted against the lattice extent L. The filled circles 
represent the lattice QCD data Eq of the lowest state in (I,J P ) = (0, j ) channel. The open symbols are 
the sum E$~ l + E^~ l of energies of nucleon E$~ l and Kaon £^"' 1 with the smallest lattice momentum 

|p| = ^|n| = 2n/L. The upper line represents ^Jm^+ |p| 2 + y^Mf~+ |p| 2 with |p| = 2k /L the smallest 
relative momentum on the lattice. The lower line represent the simple sum +Mk of the masses of 
nucleon Mn and Kaon Mk- We adopt the central values of Mn and Mk obtained on the largest lattice to draw 
the two lines. 

The filled circles in Fig. [I] show the lowest-state energies Eq in / = and J P = \ channel on 
four different volumes. Here the horizontal axis denotes the lattice extent L in the lattice unit and 
the vertical axis represents the energy of the state. The lower line denotes the simple sum M# +M& 
of the nucleon mass and Kaon mass Mk obtained with the largest lattice. We here simply use 
M N + M K (L = 16) as a guideline. 

As for the lowest state of this channel, we find that the energy of the state takes almost constant 
value against the volume variation and coincides with the simple sum M^+Mk- It is then concluded 
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that the lowest state in / = and J p = A channel is the NK scattering state with the relative 
momentum |p| =0. The good agreement with the sum M^ + Mk implies the weakness of the 
interaction between N and K. 

The (I,J P ) = (0,±~) state is one of the candidates for ©+(1540). Since ©+(1540) is located 
above the NK threshold, it would appear as an excited state in this channel. The volume dependence 
of the energy of each state can be used to distinguish a possible resonance state from NK scattering 
states; it is expected that the energies of resonance states have small volume dependence, whereas 
the energies of NK scattering states are expected to scale according to the lattice size L. 

A possible candidate for the volume dependence of the energies of NK scattering states is the 
simple formula as E" K (L) = ^M^ + |^n| 2 + \Jm^+ |^n| 2 with the relative momentum ^n be- 
tween N and K in finite periodic lattices. There may be some corrections to £'™ K (L) in practice. We 
therefore estimate three possible corrections [Q] ; the existence of the NK interaction, the application 
of the momenta on a finite discretized lattice and the estimation of the implicit finite-size effects. 
We find that these corrections lead to at most a few % deviations from E^ K (L) and we neglect these 
corrections for simplicity and use the simple form E" K (L) in the following discussion. 

We now compare the lattice data with the expected behaviors -E^ -1 for the 2nd-lowest 
NK scattering states. The filled squares in Fig. [j] denote E^, the 2nd-lowest-state energies in this 
channel. (The black and gray symbols are the lattice data obtained by the fits with two different fit 
ranges, which help the readers to see the fit-range dependences.) The upper line shows E^ -1 es- 
timated with the next-smallest relative momentum between N and K, and with the masses M K and 
Mn extracted on the L = 16 lattices. Although the lattice QCD data £*{" and E'nk" 1 almost coincide 
with each other on the L = 16 lattices, which one may take as the characteristics of the 2nd-lowest 
scattering state, the data E^ do not follow E^ -1 in the smaller lattices. Especially when the quarks 
are heavy, composite particles would be compact objects and show smaller finite volume effects be- 
sides those arising from the lattice momenta p = ^n. The statistical errors are also well controlled 
for the heavy quarks. Thus, the significant deviations in 1.5 < L < 3 fm with the combination of the 
heavy quarks, such as (k u j, k s ) = (0.1600,0.1600), are reliable and the obtained 2nd-lowest states 
are difficult to be explained as the NK scattering states. One can understand this behavior with 
the view that the 2nd-lowest state is a resonance state rather than a scattering state. In fact, while 
the data with the lighter quarks have rather strong volume dependences which can be considered 
to arise because of the finite size of a resonance state, the lattice data exhibit almost no volume 
dependence with the combination of the heavy quarks especially in 1.5 < L < 3 fm, which can be 
regarded as the characteristic of resonance states. 

3.2 volume dependence of weightfactors 

For further confirmation, we investigate the volume dependence of the spectral weight [j|]. 
In the case when a correlation function is constructed from a point-source and a zero-momentum 
point-sink, the weight factor Wi takes an almost constant value if |/) is the resonance state where the 
wave function is localized. If the state |/) is a two-particle state, the situation is more complicated. 
Nevertheless, when the interaction between the two particles is weak, which is the case in the 
present situation, the weight factor is proportional to y. When a source is a wall-type operator 
®waii(fsrc) as taken in this work, a definite volume dependence of Wj is not known. Therefore, we 



071/4 



Search for the S=+l pentaquarks in quenched lattice QCD 



Tom T. Takahashi 



re-examine the lowest state and the 2nd-lowest state in (I,J P ) = (0, \ ) channel with the hopping 
parameter (k u j,k s ) = (0.1600,0.1600) and the locally-smeared source ® 2 sm ear(fsrc) [[§], which is 
introduced to partially enhance the ground-state overlap. We additionally employ 14 3 x 24 lattice 
for this analysis. We extract the weight factors Wo and W\ using the two-exponential fit as 



0.0015 



g 0.001 

% 0.0005 
3 





=(0.1600,0.1600) 




o 16-18 




□ 17-19 


I 


o 18-20 




1 I 



1000 2000 3000 4000 5000 
Volume V [lattice unit] 



0.004 

g 0.003 
o 

3 0.002 
| 0.001 






(«i,j,Kj=(0. 1600,0. 1600) 




! ^ a 16-18 




< 17-19 




t v 18-20 



1000 2000 3000 4000 5000 
Volume V [lattice unit] 



Figure 2: The left figure shows Wo for the lowest state in (I,J P ) = (0, \ ) channel and the right figure 
shows W\ for the 2nd-lowest state. In the case when the weight factor W,- for the ;-th state |/) in a point-point 
correlator shows no volume dependence, ;) is likely to be a resonance state. On the contrary, when the 2-th 
state |/) is a two-particle state, W; scales according to 1 /V. 



I 5 (@ 2 (x,r + f src )@ 2 smear (0,f src )} = W e' E o T +Wye~ E i T , with the exponents E Q and£f fixed. We 
try the various fit ranges as (r m j n ,r max )= (16,18), (17, 19), (18,20) to see the fit-range dependences. 
Fig. ^includes all the results and one can find that the global behaviors are almost the same among 
these three. 

The figures in Fig. || show the weight factors Wq and W\ of the lowest and next-lowest state in 
(/, J p ) = (0, \ ) channel, respectively. We find that the volume dependence of Wo on V is consistent 
with 77 , which is expected in the case of two-particle states. On the other hand, W\ shows almost 
no volume dependence against V, which is the characteristic of the state in which the relative wave 
function is localized. This result can be considered as another evidence of a resonance state lying 
slightly above the NK threshold. 



4. Lattice QCD Result in (IJ P ) = (0, \ ) channel 

In the same way as (I,J P ) = (0,i ) channel, we have attempted to diagonalize the correla- 
tion matrix in (I,J P ) = (0, ^ + ) channel using the wall-sources & wa n(t) and the zero-momentum 
point-sinks Lx®(^,0- In this channel, the diagonalization is rather unstable and we find only one 
state except for tiny contributions of possible other states The data have almost no volume 
dependence and they coincide with the solid line which represents the simple sum M^* + Mk of 
Mjy and Mk, with M^* the mass of the ground state of the negative-parity nucleon. The state we 
observe is therefore concluded to be the N*-K scattering state with the relative momentum |p| =0. 



5. Summary 

We have studied the (S,I,J) = (+1,0, \) states on 8 3 x 24, 10 3 x 24, 12 3 x 24 and 16 3 x 24 
lattices at /3=5.7 using quenched lattice QCD. We have performed a systematic study by (1) the 
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variational analysis using correlation matrices with very high statistics and (2) the study of the 
volume dependences of eigenenergies and spectral weight factors with several different physical 
volumes in a proper way for the first time. From the conelation matrix of the operators, we have 
successfully obtained the energies of the lowest state and the 2nd-lowest state in the (I,J P ) = 
(0, - ) channel. The volume dependence of the energies and spectral weight factors show that 
the 2nd-lowest state in this channel is likely to be a resonance state located slightly above the NK 
threshold and that the lowest state is the NK scattering state with the relative momentum |p| = 0. 
As for the (I,J P ) = (0, 5 ) channel, we have observed only one state in the present analysis, which 
is likely to be a N*K scattering state of the ground state of the negative-parity nucleon N* and Kaon 
with the relative momentum |p| =0. 

The chiral-extrapolated value of the eigenenergy is j~ of the 2nd-lowest state in the negative 
parity channel is 1.500(52) in the lattice unit, which is about 120% of the NK threshold 1.227(6). 
This state is rather heavier than the mass of + (154O) in the real world. One possibility is the 
systematic errors from the discretization, the chiral extrapolation, or quenching. Another possibility 
is that the observed 2nd-lowest state might be a signal of a resonance state lying higher than + in 
the quenched QCD. It is difficult to give a clear explanation at this point and more extensive studies 
on finer lattices with lighter quark masses in unquenched QCD will be required. However, we can 
at least conclude that our quenched lattice calculations suggest the existence of a resonance-like 
state slightly above the NK threshold for the parameter region we have investigated. 
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